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ON  THE  PROSPECTS  FOR  ARTIFICIALLY  INDUCING  EQUATORIAL  SPREAD  F 


1.  Introduction 

The  Brazil  Ionospheric  Modification  Experiment  [BIME)  is  a  project 
sponsored  by  the  Air  Force  Geophysics  Laooratory  [Narcisi,  1982],  the  aim  of  \ 

which  is  to  produce  a  perturbation  in  the  electron  density  of  the  equatorial 
ionosphere  of  sufficient  magnitude  to  trigger  (i.e.,  provide  a  seed  for)  the 
collisional  Rayleigh-Taylor/ExB_  gradient  drift  instability,  the  instability 
believed  to  be  responsible  (see  the  theory  review  article  by  Qssakow,  1981  and 
references  therein)  for  naturally  occurring  equatorial  spread  F  (ESF).  The 
aim  of  this  experiment  is  to  make  possible  the  observation  of  the  evolution  of 
this  instability  under  more  controlled  circumstances  than  those  characterizing 
naturally-occurring  ESF.  In  this  report  we  shall  attempt  to  address  each  of 
the  factors  which  might  bear  on  the  success  or  failure  of  such  an  experiment, 
and  to  suggest  ways  by  which  the  probability  of  experiment  success  might  be 
maximized. 

In  Figure  1  we  show  the  geometry  of  the  physical  problem  of  interest  (see 
also  Zalesak  et  al. ,  1982).  The  BIME  experiment  involves  a  rocket -launched 
chemical  release  in  the  equatorial  F  region  ionosphere,  the  design  of  which  is 
to  deplete  or  enhance  the  electron  density  around  the  release  point.  Remote 
sensing  and  other  rocket  probes  would  then  monitor  the  progress  of  the 

(hopefully)  induced  instability. 

i 

\ 

2.  Theory 

For  the  analysis  in  this  section,  and  for  the  numerical  simulations  to  be 

presented  later,  we  make  the  assumption  that  the  physical  state  depicted  in 

Figure  1  can  be  accurately  modeled  by  straightening  the  magnetic  field  lines 
and  by  representing  the  distribution  of  plasma  along  magnetic  field  lines  as 
an  array  of  three  distinct  thin  layers  of  plasma  connected  by  magnetic  field 
lines,  as  depicted  in  Figure  2.  The  center  layer  represents  the  equatorial 
nighttime  F  region  plasma,  while  the  upper  and  lower  layers  represent  the 
remaining  northern  and  southern  hemisphere  plasma  respectively,  including  the 
E  region  plasma.  This  model  is  described  in  detail  in  Zalesak  et  al. 

[1982].  At  this  juncture  it  should  be  noted  that  our  previous  study  on 
artificially  created  equatorial  spread  F  (Qssakow  et  al.  ,  1978)  considered 
only  the  center  layer.  Thus,  the  effects  of  plasma  at  lower  altitudes  away 
from  the  the  equatorial  plane,  but  on  magnetic  field  lines  threading  the 
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equatorial  region,  on  the  bubble  dynamics  was  not  considered  (in  reality  that 

model  was  valid  only  when  the  depletion  and  not  the  background  ionospheric 

integrated  density  and  Pedersen  conductivity  dominate  the  dynamics).  We  now 
list  here  the  primary  assumptions  in  the  present  model:  1)  the  electric 

fields  JE  of  interest  are  electrostatic  in  nature  and  hence  derivable  from  a 
scalar  potential:  2)  the  conductivity  along  magnetic  field  lines 

Is  sufficiently  high  that  we  may  treat  the  magnetic  field  lines  as  equipoten- 
tials:  $  ■  $(x,y)  (see  Figure  2);  3)  currents  between  layers  are  carried 

solely  by  electrons,  i.e. ,  there  is  no  ion  tranport  along  magnetic  field 
lines;  4)  layers  1  and  3  (see  Figure  2)  are  initially  uniform  in  both 

electron  density  and  Pedersen  conductivity,  and  remain  so  during  all  times  of 
interest,  enabling  us  to  treat  them  as  a  passive  electrically  conducting  load 
on  our  system.  For  further  details,  the  reader  is  referred  to  Zalesak  et  al. 
[1982]. 

Within  the  context  of  the  above  model,  the  JjxB_  gradient 
drif t/collisionally  dominated  Rayleigh-Taylor  instability  can  be  triggered 
whenever  a  certain  geometric  relationship  holds  between  the  electron  density 
gradient  in  the  equatorial  ionosphere  and  the  forces  acting  on  the  plasma.  In 
particular,  an  infinitesimal  perturbation  on  the  system  will  grow 

Y  t 

exponentially,  i.e.,  ap  e  where  t  is  time,  with  growth  rate  y  given  by 
(Zalesak  et  al.  ,  1982) 
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where  and  B_  are  the  electric  and  magnetic  fields  respectively,  is  the 

neutral  wind  velocity,  a  is  the  gravitational  acceleration,  Z  .  is  the 

pi 

magnetic  field  line  integrated  Pedersen  conductivity  of  layer  i(i  -  1,2,3),  n2 
is  the  electron  density  in  layer  2,  and  is  the  ion-neutral  collision 
frequency  in  layer  2.  R  is  a  term  due  to  recombination  chemistry,  which  will 
be  discussed  shortly.  In  Figure  2  it  can  be  seen  that  Vn2  for  our  unperturbed 
ionosphere  is  in  the  y  direction  (vertical  in  the  equatorial  plane).  Since 
the  vertical  component  of  the  neutral  wind  is  in  general  quite  small  with 
respect  to  the  other  terms  in  the  growth  rate,  we  may  neglect  this  term  in 
Eq.  (1). 


The  presence  or  non-presence  of  the  term  R  in  Eq.  (1)  is  a  function  of 
ones  choice  of  a  zeroth  order  profile  for  the  equatorial  ionosphere*  In  the 
original  papers  of  Scannapieco  and  Ossakow  [1976]  and  Ossakow  et  al*  [1979], 
the  continuity  equation  for  the  electrons  was  of  the  form 

3n2/3t  +  V*(n2v)  *  -  v  (n2-n  )  (2) 

where  ^  is  the  electron  velocity.  The  recombination  coefficient  represents 
the  sum  of  the  rate-limiting  charge  exchange  reactions  of  0+  with  molecular 
oxygen  and  nitrogen.  In  this  equation  nQis  the  assumed  equilibrium  electron 
density  in  the  equatorial  ionosphere,  i*e* ,  8nQ/3t  *  0.  Since  recombination 
chemistry  itself  can  only  result  in  a  loss  of  electrons  at  a  rate  -v^n2,  we 
note  that  Eq.  (2)  implicitly  assumes  the  presence  of  a  source  of  ionization  in 
the  equatorial  ionosphere  of  value  +  vRn0*  The  presence  of  such  a  true 
ionization  source  in  the  nighttime  equatorial  ionosphere  is  subject  to  some 
question  and  subsequent  work  by  Zalesak  et  al.  [1982]  dropped  the  assumption 
of  the  existence  of  such  a  term,  resulting  in  an  electron  density  equation  of 
the  form 

3n2/3t  +  V«(n£v)  ■  -  vRn2  (3) 

Note  that  Eq.  (3)  does  not  yield  a  time-invariant  equilibrium  for  the 
unperturbed  electron  de'vity  profile,  but  rather  a  profile  which  steadily 
decays  through  the  night. 

A  linear  stability  analysis  of  the  full  set  of  equations  describing  the 
dynamics  of  the  ExB_  gradient  drif t/collisionally  dominated  Rayleigh-Taylor 
instability  in  the  equatorial  ionosphere  [see  Ossakow  et  al.  ,  1979],  shows 

that 

v  for  Eq.  (2) 

R  -  {  (4) 

0  for  Eq.  (3) 

Thus  the  effect  of  recombination  chemistry  on  the  instability  growth  rate  is 
determined  by  ones  choice  of  a  model  for  :he  ambient  ionosphere. 
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Equation  (1),  together  with  a  choice  of  the  ambient  ionospheric  state,  as 
just  discussed,  provides  us  with  a  criterion  for  determining  whether  or  not  an 
infinitesimal  perturbation  to  the  ionosphere  will  be  amplified;  but  since  the 
experiment  is  not  constrained  to  infinitesimal  perturbations,  one  is  led  to 
consider  both  of  the  following  scenarios: 

A)  the  ambient  ionosphere  is  stable  to  infinitesimal  perturbations 

B)  the  ambient  ionosphere  is  unstable  to  infinitesimal  perturbations 

Case  A,  that  of  a  linearly  stable  ionosphere,  implies  that  ones  only  hope 
of  producing  artificial  ESF  is  to  perturb  the  ionosphece  with  sufficient 
amplitude  to  drive  the  dynamics  fully  into  the  nonlinear  regime  immediately. 
A  series  of  numerical  simulations  could  be  used  to  analyze  this  case 
thoroughly,  but  the  results  will  undoubtedly  depend  strongly  on  both  the 
degree  of  stability  and  on  the  mechanism  providing  the  stabilization,  as  well 
as  on  how  hard  the  ionosphere  is  driven  by  the  chemical  release. 

Case  B,  that  of  a  linearly  unstable  ionosphere,  immediately  leads  one  to 
ask:  why  has  ESF  not  already  occurred?  In  fact,  one  Is  generally  tempted  Co 
treat  the  non-occurrence  of  ESF  as  an  indication  of  the  stability  of  the 
ionosphere.  In  this  case  one  is  hoping  to  find  an  ionosphere  which  has  only 
been  subject  to  extremely  small  amplitude  perturbations,  or  one  whose  growth 
rate  is  positive  but  relatively  small,  or  both,  in  which  instance  he  can  be 
reasonably  sure  of  triggering  the  instability  with  the  chemical  release. 
However,  it  may  be  difficult  to  prove  a  posteriori  that  ESF  would  not  have 
occurred  even  in  the  absence  of  the  release. 

In  either  case  A  or  B  above,  one  would  be  advised  to  monitor  the 
ionosphere  using  ionosondes  and  other  instrumentation  for  many  nights  prior  to 
the  launch  in  order  to  get  a  "feel"  for  the  statistics  of  natural  ESF 
occurrence* 

It  is  our  view  that  case  A  above  represents  an  untenable  approach  to  this 
experiment.  In  addition  to  the  problems  already  mentioned,  the  chances  of 
driving  the  dynamics  fuly  into  the  nonlinear  regime  via  the  chemical  release 
would  appear  to  be  quite  slim  with  existing  payloads  (such  as  those  in  BIME, 
considered  by  Narcisl,  1982).  The  physics  of  the  instability  depends  on  the 
total  magnetic  field  line  Integrated  Pedersen  conductivity  as  well  as  the 
field  line  integrated  electron  density  in  level  2  (not  simply  the  electron 
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density  exactly  at  the  magnetic  equator)*  Anderson  and  Bernhardt  [1978]  have 
shown  that  even  a  10  kg  H2  release  (a  large  release)  at  350  km  altitude  will 
yield  only  a  maximum  of  6%  reduction  in  total  field  line  integrated  Pedersen 
conductivity.  The  task  of  finding  an  unstable  but  as  yet  unstructured 

ionosphere  (case  B  above)  would  appear  to  be  much  easier.  However  (and  we 
cannot  stress  this  point  too  strongly)  one  must  be  certain  that  he  has  not 
misidentif led  a  stable  ionosphere. 

3.  Numerical  Simulations 

Details  of  the  numerical  techniques  used  in  this  simulation  and  of  the 
general  computational  procedure  may  be  found  in  Qssakow  et  al.  [1979]  and 
Zalesak  et  al.  [1982]  and  in  the  references  therein.  Briefly,  the  numerical 
calculations  to  be  presented  were  performed  on  a  two-dimensional  Cartesian 

mesh  using  40  points  in  the  x  (east-west)  direction  and  140  points  in  the  y 
(vertical  at  the  equator)  direction.  The  uniform  grid  spacing  was  3  km  in  the 
y  direction  and  5  km  in  the  x  direction.  Periodic  boundary  conditions  were 
imposed  on  n,  n2  and  on  the  electrostatic  potential  <p  in  the  x  direction.  At 
the  top  boundary  in  y  we  set  3n/3y  =  3n2/3y  *  3<j>/3y  *  0,  while  at  the  bottom 
boundary  we  set  3n/3y  -  3n2/3y  =*  <{>  »  0.  Here  n  is  the  electron  density 

exactly  at  the  magnetic  equator,  to  be  distinguished  (see  below)  from  n2,  the 
average  electron  density  in  layer  2.  In  principle,  each  of  the  three  layers 
must  be  represented  on  a  separate  2-D  mesh,  but  this  need  not  be  done  here 
since  one  of  our  assumptions  Is  that  layers  1  and  3  are  uniform  and  remain  so 
during  the  course  of  the  calculation  and  since  <(>  is  identical  on  each  mesh. 
Both  the  continuity  equation  (of  the  form  of  Eq.  3  rather  than  Eq.  2)  and  an 
elliptic  equation  for  the  electrostatic  potential  <f>  are  solved  on  this  grid 
(see  Zalesak  et  al.  a  1982).  Realistic  values  for  the  ion-neutral  collision 
frequency  and  recombination  coefficient  vR  in  layer  2  were  used,  and  are 
given  in  Qssakow  et  al.  [1979].  Layers  1  and  3  were  each  assigned  an 

integrated  Pedersen  conductivity  equal  to  6%  of  the  maximum  of  that  found  in 
layer  2  for  our  "unstable”  case  B  (to  be  described  shortly). 

The  form  of  the  ionospheric  hole  caused  by  the  chemical  release  was  taken 

to  be 


n(x,y)/no(y)  -  1  -  0.95  exp[ (r/37)"2,35] 


(5) 
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where  n  is  the  electron  density  exactly  at  the  magnetic  equator,  nQ  Is  the 
ambient  electron  density  at  the  equator  in  the  absence  of  the  chemical 
release,  and  r  is  the  distance  in  kilometers  from  the  point  of  release.  Note 
that  we  distinguish  between  n,  the  electron  density  exactly  at  the  equator, 
and  n2,  the  average  electron  density  in  layer  2.  In  the  absence  of  the 
perturbation,  we  take  n  «  nj  at  t  *  0.  Hence  nQ  represents  the  ambient 
profile  for  both  n  and  n£.  The  form  and  coefficients  of  Eq.  (5)  were  chosen 
to  approximate  the  Natal  1900  LT  release  in  Fig.  5  of  Mendillo  [1981].  The 
form  of  nQ(y)  is  identical  to  that  used  in  our  previous  work  [Ossakow  et  al.  , 
1979].  Note  that  Eq.  (5)  gives  a  large  (95%)  depletion  in  electron  density  at 
the  center  of  the  hole.  However,  this  figure  is  expected  to  hold  only  exactly 
at  the  release  point  (presumably  the  magnetic  equator).  Recalling  that  our 
layer  2  is  meant  to  comprise  all  of  the  plasma  within  several  degrees  of  the 
magnetic  equator,  it  is  appropriate  to  substantially  reduce  the  effect  of  the 
release  in  terms  of  its  effect  on  the  total  magnetic  field  line  integrated 
electron  density  in  layer  2.  We  therefore  take 

^2(x,y)/no(y)  =■  1  -  0.05  exp[ (r/37)~2*35 ]  (6) 

where  ^  is  the  field-line  averaged  electron  density  in  layer  2  and  r  is  now 
two-dimensional.  The  5%  maximum  depletion  given  in  Eq.  (6)  was  chosen  to  be 
in  the  general  range  of  estimates  of  H2  releases  larger  than  BIME  given  by 
Anderson  and  Bernhardt  [1978],  and  from  Eq.  (5)  corresponds  to  layer  2  being 
approximately  1000  km  thick  (along  B_). 

In  order  to  clarify  the  roles  of  n2  and  n,  and  for  the  sake  of 
completeness,  we  enumerate  here  the  equations  which  we  solve  numerically: 


3n2/3t  +  7i*(n2l^  *  "  uRn2 

(7) 

3n/3t  +  V  |  • (nv )  ■  -  v  n 

1  *— ■  R 

(8) 

VHvinn2  +  Cb)7l*]  "  h  ^Eox  ‘  T^vinn2 

in 

(9) 

L-f-  (So  -  v>  *  2. 

(10) 
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(13) 

(14) 


Note  that  we  have  written  Eq.  (9)  to  show  explicitly  the  lack  of  any 
dependence  on  n. 

In  light  of  the  above  discussion,  plots  of  both  n  and  ri2  will  be  shown, 
with  the  understanding  that  while  it  is  n  which  is  likely  to  be  measured  with 
rocket  probes  or  other  sensors  at  the  equator,  it  is  ^  which  drives  the 

[  physics*  In  terms  of  the  numerical  calculations  to  be  shown,  n  can  be  thought 

I, 

of  as  a  tracer  fluid  which  merely  describes  what  is  happening  exactly  at  the 
equator  (neglecting  parallel  to  B_  diffusion  effects),  but  whose  effect  on  the 
dynamics  can  only  be  felt  by  integrating  along  an  effect  which  is  already 
folded  into  the  computation  of  vi2»  Two  kinds  of  plots  will  be  presented: 
contours  of  constant  n2(x,y,t);  anc*  contours  of  constant  n(x,y,t)* 
Superimposed  on  each  contour  plot  is  a  dashed  line  depicting  nQ(y)  for 
reference  purposes.  Our  nQ(y)  profile  is  such  that  the  minimum  electron 
density  scale  length  L  »  n0(3nQ/3y)  ^  is  10  km,  and  varies  in  the  runs  to  be 
presented  only  in  terms  of  the  altitude  of  the  F2  peak,  as  can  be  seen  in  the 
plots  to  be  shown* 

We  chose  to  perform  two  calculations  which  represent  the  two  possible 
scenarios  previously  discussed:  the  stable  ionosphere  (case  A)  and  the 
unstable  ionosphere  (case  B)*  Since  we  have  chosen  Eq.  (3)  as  our  continuity 
equation,  we  have  implicitly  assumed  that  recombination  chemistry  cannot 
affect  the  growth  rate  (see  Eq.  4),  so  the  only  truly  stabilizing  mechanism  we 
have  at  our  disposal  is  westward  ambient  electric  field  J^.  We  can  also 
mitigate  the  effect  of  the  positive  g/vin  term  in  (1)  by  lowering  the 
ionosphere  (which  increases  v^).  Accordingly,  our  simulations  were  run  with 
the  following  parameters 
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F2  peak  altitude 
chemical  release  altitude 
^  x  velocity 


Case  A  Case  B 
375  km  430  km 
320  km  375  km 
70  m/sec  downward  0 


Case  A  is  such  that  the  entire  ionosphere  is  linearly  stable.  Case  3  is 
such  that  large  regions  of  the  bottomside  F  region  are  linearly  unstable,  and 
in  fact  this  run  is  identical  to  calculation  2LE  in  Zalesak  et  al.  [1982] 
except  for  the  form  of  the  perturbation.  Figures  3-6  show  the  results  of  the 
simulations  at  selected  times.  Plots  of  both  n  and  n£  are  shown.  Note  in 
case  A  that  although  the  ionosphere  is  changing  (due  to  the  downward  ^  x  B 
velocity  and  recombination  chemistry)  the  perturbation  is  in  fact  de  ^ed.  Also 
note  that  the  perturbation  in  n  will  never  totally  damp  even  though  at  in  n2 
(the  one  driving  the  dynamics)  will.  In  case  E  we  see  that  wa  indeed 

initiated  artificial  ESF  with  the  chemical  release.  The  importa*  hing  to 
note  is  that  linear  stability  theory  is  an  accurate  predict  the 

experimental  success  or  failure,  especially  when  the  chemical  release  cannot 
reach  the  nonlinear  regime  immediately. 

4.  Finding  an  Unstable  Ionosphere 

The  primary  finding  of  the  results  of  the  numerical  simulations  presented 
in  Section  3  is  that  they  bear  out  the  validity  of  a  linear  stability  analysis 
in  terms  of  analyzing  the  BIME  experiment.  Accordingly,  we  conclude  that  the 
most  viable  approach  to  take  in  the  BIME  experiment  is  to  find  an  unstable, 
but  as  yet  unstructured  ionosphere.  As  mentioned  before,  the  drawback  to  this 
approach  is  the  possibility  of  having  misidentif ied  a  stable  ionosphere. 
Therefore  let  us  appeal  to  Eq.  (1)  and  to  other  known  stabilizing  influences 
in  the  equatorial  ionosphere  and  draw  up  a  list  of  possible  stabilization 
mechanisms  of  which  we  must  beware: 

1)  recombination  chemistry,  if  we  believe  there  to  be  a  true  source  of 
ionization  present  (Eq.  (2))  such  as  to  keep  nQ(y,t)  «  nQ(y,t*0); 

2)  a  westward  ambient  electric  field; 

3)  shear  stabilization,  caused  by  a  vertical  shear  in  the  horizontal 
plasma  motion  of  the  ambient  ionosphere  [Perkins  and  Doles,  1975] ;and 
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4)  inordinately  small  positive  growth  terms  in  Eq.  (1),  such  as  large 

small  gradients  in  n2»  or  a  large  background  conductivity 

L  .  +  Z 
Pi  p3 

Optimization  of  31ME  success  probability  can  be  reduced  to  minimizing  the 

impact  of  each  of  the  above  mechanisms.  If  we  believe  recombination  chemistry 

to  be  a  problem,  then  since  we  know  that  v  decreases  with  altitude,  we 

conclude  that  we  should  only  launch  Into  ionospheres  which  are  reasonably  high 

in  altitude.  A  westward  electric  field  can  be  avoided  by  noting  that  such  a 

field  would  produce  downward  plasma  velocities.  An  ionosonde  with  doppler 

capabilities  could  detect  such  an  occurrence  and  generate  a  "no-go"  signal. 

The  shear  stabilization  mechanism  is  not  well  understood,  but  it  should  be 

noted  that  large  plasma  shears  are  an  indication  of  large  background 

conductivities  [Zalesak  et  al. ,  1982],  and  hence  should  be  avoided  on  both 

counts.  Large  values  of  v  can  be  avoided  by  the  same  mechanism  as  that  of 

avoiding  recombination  chemistry,  i.e.  ,  launching  only  into  reasonably  high 

ionospheres,  since  v,  ,  like  v  .  decreases  with  altitude.  Ionosondes  can  be 

in  R 

used  to  measure  gradients  in  n,  which  presumably  reflect  gradients  in  for 
the  undisturbed  ionosphere.  As  was  already  mentioned,  although  background 
conductivities  cannot  be  measured  directly,  the  degree  of  plasma  shear  can  be 
used  to  give  an  indication  of  their  size.  Better  yet  would  be  a  measure  of 
the  difference  between  the  plasma  and  neutral  gas  zonal  velocity  (see  Eq.  38 
in  Zalesak  et  al.  [1982]). 

5.  Conclusions  and  Recommendations 

It  is  our  conclusion  that  the  success  of  the  B IME  experiment  (and  any 
other  that  cannot  reach  the  nonlinear  regime  immediately)  depends  crucially  on 
one's  ability  to  find  an  unstable  but  as  yet  unstructured  ionosphere.  This  is 
not  a  trivial  task.  The  only  means  of  truly  ensuring  that  this  is  the  case 
involves  in-situ  measurements  of  several  physical  quantities  all  along 
magnetic  field  lines  and  throughout  the  equatorial  ionosphere.  Realistically, 
the  best  one  can  do  is  to  try  to  minimize  the  risk  by  avoiding  those  specific 
situations  which  would  indicate  ionospheric  stability.  Specifically,  we 
recommend  the  following: 

1)  Try  to  launch  into  a  high  ionosphere,  In  order  to  maximize  the 
terra  g/v^  In  the  linear  growth  rate,  and  to  minimize  the  stabilizing 
effects,  if  any,  of  recombination  chemistry. 
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2)  Try  to  launch  into  an  ionosphere  that  is  rising,  or  at  the  very  least 

one  whi  i  is  not  falling,  in  order  to  minimize  the  chance  of 

encountering  a  stabilizing  ambient  westward  electric  field.  An 
ionosonde  with  doppler  capabilities  would  be  very  useful  for  this 
purpose. 

3)  Try  to  avoid  launching  into  regions  of  large  horizontal  shear  in 

plasma  motion. 

4)  Try  to  have  available  in  the  field  the  ability  to  measure  logarithmic 

electron  density  gradients,  compute  the  approximate  altitude  of 

maximum  growth  from  Eq.  (1),  and  perform  the  chemical  release  of  that 
altitude. 

5)  Monitor  the  ionosphere  using  ionosondes  for  many  nights  prior  to 

launch  to  get  a  "feel"  for  the  statistics  of  natural  equatorial 
spread  F  occurrence. 


I 


EQUATORIAL 
SPREAD  F  MODEL 


4000 

km 


1125  km 


Fig.  1  —  Diagram  of  the  equatorial  ionosphere  and  of  the  neighboring 
regions  which  have  physical  relevance  to  equatorial  spread  F  (ESF) 
processes,  including  the  E  region  plasma  at  higher  and  lower  latitudes. 
These  regions  are  electrically  coupled  to  the  equatorial  F  region  iono¬ 
sphere  by  the  high  conductivity  along  magnetic  field  lines.  Plasma  is 
actually  distributed  all  along  these  field  lines,  but  in  this  study  we 
shall  make  the  assumption  that  this  system  can  be  modeled  accurately 
by  three  planes  of  plasma  connected  by  straight  field  lines,  as  shown 
in  Fig.  2.  One  of  these  three  layers  (layer  2  in  Fig.  2)  is  shown  here 
as  the  “computational  plane.” 
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Fig.  2  —  The  “three  layer”  model  of  the  physical  system  depicted  in  Fig.  1. 
All  plasma  in  the  vicinity  of  the  equatorial  plane  has  been  compressed  into 
layer  2,  while  the  remaining  northern  and  southern  hemisphere  plasma  has 
been  compressed  into  layers  1  and  3  respectively.  Further,  the  magnetic  field 
lines  have  been  straightened  so  we  can  deal  in  cartesian  coordinates  x,  y,  and 
z  as  shown  in  the  figure.  The  plasma  in  layers  1  and  3  is  assumed  to  be  uni¬ 
form  and  free  of  any  external  driving  force.  The  equatorial  layer  2  is  assigned 
a  realistic  initial  distribution  of  electron  density  N0(y),  and  ion-neutral  colli¬ 
sion  frequency.  In  addition,  gravity  points  in  the  negative  y  direction. 
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Fig.  3  Contours  of  constant  n  for  case  A  (stable  ionosphere)  at  various  times  in 
the  calculation.  At  time  0  seconds  the  contours  are  labeled  in  units  of  electrons/ 
cm3  in  E  format  notation  (1.0  El  =  1  x  10*.  etc.).  The  unperturbed  ionosphere 
was  initially  laminar  (independent  of  x,  the  east-west  direction)  and  is  exhibited  by 
the  dashed  curve  showing  nQ(y).  This  curve  is  labeled  at  the  top  of  the  figure.  The 
perturbation  shown  at  0  seconds  is  described  by  Eq.  (5).  Two  full  periods  of  the 
calculation  are  shown  to  emphasize  the  periodic  boundary  conditions  used  in  the 
x-direction.  The  observer  is  looking  southward  so  that  B  is  out  of  the  figure. 
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OICY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
OICY  ATTN  NRT 

OICY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MX 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

OICY  ATTN  MNNL  LTC  KENNEDY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EF.P  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 

OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 

OICY  ATTN  DOC  CON  FOR  E.  JONES 

OICY  ATTN  DOC  CON  FOR  J.  MALIK 

OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

OICY  ATTN  DOC  CON  FOR  J.  ZINN 

OICY  ATTN  DOC  CON  FOR  P.  KEATON 

OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  1HORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 

OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

(ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

OICY  ATTN  R.  MOORE 
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INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

OlCY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

OlCY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
OlCY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OlCY  ATTN  R.  GRUBB 
OlCY  ATTN  AERONOMY  LAB  G.  REID 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OlCY  ATTN  I.  GARFUNKEL 
OlCY  ATTN  T.  SALMI 
OlCY  ATTN  V.  JOSEPHSON 
OlCY  ATTN  S.  BOWER 
OlCY  ATTN  N.  STOCKWELL 
OlCY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OlCY  ATTN  RADIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.O.  BOX  983 
BERKELEY,  CA  94701 
OlCY  ATTN  J.  WORKMAN 
OlCY  ATTN  C.  PRETTIE 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

OlCY  ATTN  G.  KEISTER 
OlCY  ATTN  D.  MURRAY 
OlCY  ATTN  G.  HALL 
OlCY  ATTN  J.  KENNEY 

BROWN  ENGINEERING  COMPANY,  INC. 

CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OlCY  ATTN  ROMEO  A.  DELIBERIS 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 

P.O.  BOX  6049 

SAN  DIEGO,  CA  92106 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OlCY  ATTN  D.B.  COX 
OlCY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  M D  20734 
OlCY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OlCY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OlCY  ATTN  H.  LOGSTON 
OlCY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OlCY  ATTN  J.  ROBERTS 
OlCY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OlCY  ATTN  M.H.  BORYNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OlCY  ATTN  F.  RE  I  BERT 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OlCY  ATTN  G.  MILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARBARA,  CA  93111 
OlCY  ATTN  JOHN  ISE,  JR. 

OlCY  ATTN  JOEL  GARBARINO 


1 


i 


I 


i 

! 

! 

i 


1 


23 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVAN LA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  MARSHALL  CROSS 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAURAL,  HD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 


01CY  ATTN  T.  MEAGHER 
KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS LAC 
01CY  ATTN  TIM  STEPHANS 
01CY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
0  ICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  &  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGKLIN 
OICY  ATTN  D.  CLARK 

KCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  S2647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 
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01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 
MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN- L.  GUTSCHE 
01CY  ATTN  D.  SAPPENFIELD 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOULE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD, 

SANTA  MONICA,  CA  90404 
01CY  ATTN  E.C.  FIELD,  JR, 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
0  ICY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


01CY 

ATTN 

DONALD  NEILSON 

01CY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

01CY 

ATTN 

L,L,  COBB 

OICY 

ATTN 

DAVID  A*  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

M.  BARON 

OICY 

ATTN 

RAY  L.  LEAD AB RAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

J*  PETERSON 

OICY 

ATTN 

R.  HAKE,  JR, 

OICY 

ATTN 

V,  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
01CY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R,  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D«  DEE 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
01CY  ATTN  J,  CARPENTER 
01CY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  COO  WAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 


PRINCETON  UNIVERSITY 

UNIVERSITY  OF  PITTSBURGH 

PLASMA  PHYSICS  LABORATORY 

PITTSBURGH,  PA  15213 

PRINCETON,  NJ  08540 

DR.  N.  ZABUSKY 

DR.  F.  PERKINS 

DR.  M.  BIONDI 

DR.  E.  OVERMAN 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 

UTAH  STATE  UNIVERSITY 

LA  JOLLA,  CA  92037 

4TH  AND  8TH  STREETS 

DR.  D.A.  HAMLIN 

LOGAN,  UTAH  84322 

DR.  L.  LINSON 

DR.  R.  HARRIS 

DR.  E.  FRIEMAN 

DR.  K.  BAKER 

DR.  R.  SCHUNK 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 


U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
M.  PONGRATZ 
D.  SIMONS 
G.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 
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